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We study the momentum distribution of the electrons in an extended periodic Anderson model, 
where the interaction, Ucf, between itinerant and localized electrons is taken into account. In the 
symmetric half-filled model, due to the increase of the interorbital interaction, the / electrons become 
more and more delocalized, while the itinerancy of conduction electrons decreases. Above a certain 
value of Ucf the / electrons become again localized together with the conduction electrons. In the 
less than half-filled case, we observe that Ucf causes strong correlations between the / electrons in 
the mixed valence regime. 
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I. INTRODUCTION 

Heavy-fermion and mixed valence systems are still ac¬ 
tive research fields in spite of the major achievements 
of the past few decadesf^i^ The discovery of a new crit¬ 
ical point in the pressure-temperature phase diagram of 
CeCu 2 Ge 2 and CeCu 2 Si 2 has attracted much attention 
both experimentally^r— and theoreticallyi^r— It is be¬ 
lieved that the appearance of the new critical point is 
due to the critical valence fluctuations of the Ce ion. The 
simplest model, which contains the essential physics of 
rare-earth compounds is the periodic Anderson model: 

UpAM = -t 
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where the notation is standard and W = At \s taken as 
the energy unit. It is known, however, that the mixed- 
valence regime appears always in this model as a smooth 
crossover, and valence fluctuations do not become critical 
for any choice of the parameters. A local Coulomb in¬ 
teraction between the conduction and localized electrons 
is needed for the appearance of a sharp transition and 
critical valence fluctuations Therefore we consider the 
following Hamiltonian: 

'H = fipAM + Ucf X ' (2) 

3,<r,(7' 

Previous studies revealed how Ucf affects the mixed va¬ 
lence regime and it has been shown that a first-order 
valence transition and a quantum critical point may ap¬ 
pear due to The effect of Ucf in the Hondo 

regime has been also addressed both in infinit o^^d^ and 
one spatial dimensions.— Namely, in infinite dimensions 
the symmetric model for small hybridization (V <C W) 
displays antiferromagnetic order for small Ucf which, 


however, disappears for large Ucf and charge order de¬ 
velops. In contrast, there is no such phase transition 
in one dimension due to the enhanced quantum fluctu¬ 
ations, however, for small and large Ucf the spin-spin 
and density-density correlation function, respectively ex¬ 
hibits the slowest decay. Between these two regimes there 
is a narrow region, where the local singlet formation is 
significantly enhanced. 

Our goal in this paper is to investigate the momen¬ 
tum distribution of the electrons in one dimension. It 
is known, that in higher dimensions they exhibit a jump 
at the Fermi momentum, whose size can be used to ex¬ 
tract the energy dependence of the self-energies, from 
which the many-body enhancement factor of the effec¬ 
tive mass can be obtained. Although in one dimension 
there is no such jump at the Fermi momentum, just 
a sharp change, they provide direct information about 
the spatial distribution of the electrons and the con¬ 
tent of conduction and /-electron states in the quasi¬ 
particle bands, while the previous quantum information 
analysis2£ gave only an indirect description of these quan¬ 
tities. We address the question how they are modihed by 
switching on Ucf both in the integer and mixed valence 
regimes. The density-matrix renormalization-group algo¬ 
rithm (DMRG)^“— is applied, which allows the accurate 
determination of ground state properties. We have used 
the dynamic block-state selection algorith m^^i^^ in which 
the threshold value of the quantum information loss, y, 
is set a priori. We have taken y = I0“®. A maximum of 
2000 block states is needed to achieve this accuracy, and 
the largest truncation error was in the order of I0“®. We 
investigated chains up to a maximum length L = 80 with 
open boundary conditions and performed 8-12 sweeps. 

II. RESULTS AT HALF FILLING 

The non-degenerate version of the periodic Anderson 
model can hold up to Umax = 4 electrons per lattice site, 
the average number of c and / electrons per site, and 

, respectively, can vary between zero and two. The 
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filling will refer to the ratio of the total electron density 
per site {n^ n^). In what follows we consider the sym¬ 

metric half-filled model, where = 1, and calculate the 
momentum distribution of conduction and / electrons 
which are defined as 

<7 

= \'^{fLfka) ■ ( 4 ) 
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Our DMRG calculation was performed in real space, 
therefore these quantities can be obtained by Fourier 
transforming the corresponding single particle density 
matrices, namely: 

jla 

jla 

where k = 27rn/i and n = —_L/2 — 1,... ,i/2. In our 
case these are symmetric functions, therefore we consider 
only the nonnegative k values. 

Before going into the details of the numerical results, 
we briefly recall the case when Uf = Ucf = 0 which is 
easily solvable. Thereby the Hamiltonian can be diago¬ 
nalized by an unitary transformation 

ai~'’ = -VkCk + Ukfk, 

( + ) ' 

^kfki 


where ^ creates a quasiparticle in the lower 

(upper) hybridized band with mixing amplitudes: 
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in our case Sk = —‘2t cos k and the Fermi momentum is at 
the boundary of the Brillouin zone since the lower band 
is completely filled. It is easily seen that the momen¬ 
tum distribution functions provide information about the 
mixing amplitudes, namely, the portion of conduction 
and / states in the hybridized band: 


open boundary condition. We checked for short systems 
that the momentum distributions of the interacting sys¬ 
tem calculated with periodic and open boundary condi¬ 
tions are in good agreement within our error margin. 



FIG. 1. Momentum distribution functions of the conduction 
and / electrons in the noninteracting case {Uf = Ucf = 0, 
VjW = 0.1, Ef = 0). The solid and dotted lines are obtained 
from Eq. © and ®, respectively. The symbols •, ■ denote 
the DMRG results for L — 50. 

It has been pointed out^SiS^ that strong Uf leads to 
localization of the / electrons, since the number of dou¬ 
bly occupied / levels is negligible and the ground state 
is a collective singlet. This is what we see in Fig. [2] 
(a), namely, n^{k) hardly depends on k in the Kondo 
regime, while the distribution of conduction electrons is 
just slightly affected hy Uf. The quantum information 



n^'ik) = vl, 
nk{k) = ul- 
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The momentum distribution in the noninteracting sys¬ 
tem is shown in Fig. [1] compared with the DMRG results. 
The small discrepancy between the two results is at¬ 
tributed to the open boundary condition used in DMRG. 
In the following we investigate how the interactions mod¬ 
ify the above results, using again the DMRG method with 


FIG. 2. Momentum distribution functions of the conduction 
(•) and / (■) electrons for L = 50 and n = 2. Panel (a), 
(b), (c) and (d) correspond to UcfjW = 0,1.5,1.7 and 3, 
furthermore Uf/W = 3, V/W — 0.1 and £/ = —Ufl2 in all 
cases. The lines are guides to the eye. 

analysis^S showed that as Ucf is switched on more and 
more doubly occupied / sites are created, therefore the 
/ electrons become less localized in real space. Finally, 
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when Ucf is large the c and / electrons tend to avoid 
each other and the sites are occupied by two c or two 
/ electrons in an alternating fashion. Now we examine 
how these features are reflected in the momentum dis¬ 
tributions. As Ucf is switched on, the wave number de¬ 
pendence of the conduction electrons become weaker and 
weaker and therefore less itinerant in real space, while the 
distribution of the / electrons becomes more and more 
dispersive as it can be seen in Fig. [2] (b) and (c). Above 
Ucf ~ Uf/2 + W/4: both distributions hardly depend on 
the wave number as it is observed in Fig. [5] (d). That is, 
the behavior of the momentum distributions agrees well 
with the results of the entropy analysis. 


III. AWAY FROM HALF FILLING 

In the previous section we considered the half-filled 
case. Now we discuss what happens when the ground 
state is metallic, and fix the electron density at n = 1.75. 
It has been shown,— that a stable mixed valence 
regime appears around = 2 — n in the presence of 
Ucf, and strong enough Ucf leads to a first-order tran¬ 
sition between the Kondo and mixed valence states as 
Sf is varied. One can observe in Fig. [31 that for strong 
Ucf a stable mixed valence regime appears indeed, and 
the change of the valence becomes sharp. In the follow- 



FIG. 3. The /-level occupancy as a function of e/ for L = 
80, Ucf = 0 (■) and Ucf/W = 4 (•), furthermore Uf/W = 
10, V/W = 0.2. The circled data points are used in the 
comparison in Fig. [4] The lines are guides to the eye. 

ing we investigate how the momentum distribution of the 
electrons change due to Uc / in the mixed valence regime. 
This is shown in Fig. |4l where e/, indicated by the cir¬ 
cles around the data points in Fig. [31 was chosen such 
that the occupancy of the / level is nearly the same in 
the two cases. One can clearly see that both rU(k) and 
n-^ (k) change drastically around the Fermi momentum 
when Ucf = 0. This is not surprising since our system 
is a Luttinger liquid, where a logarithmic singularity is 
expected to occur at the Fermi momentum. For a finite 
Ucf a significant amount of the conduction and / elec¬ 
trons is scattered above the Fermi momentum and the 
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FIG. 4. Momentum distribution functions of the conduction 
(•) and / (■) electrons for L = 80 and n = 1.75. Panel (a), 
(b) correspond to Ucf/W = 0, effW = 0.5 and Ucf jW = 
4, EffW = —0.75, respectively, furthermore Uf/W = 10, 
VjW = 0.2 in all cases. The lines are guides to the eye. 

/c-dependence of the electron densities is significantly re¬ 
duced in both cases, which indicates a strongly correlated 
mixed valence state. It is worth noting that the distri¬ 
bution functions do not tend to zero above the Fermi 
momentum for strong Ucf. These results agree well with 
what has been obtained in infinite dimensions using the 
Gutzwiller wave functionals the main difference is the ab¬ 
sence of the discontinuity at the Fermi momentum due 
to the one-dimensional property of the model. 


IV. CONCLUSIONS 

We have investigated an extended periodic Anderson 
model with an additional Coulomb interaction using the 
DMRG algorithm to better understand its effect on the 
momentum distribution of the electrons. In the half- 
filled, symmetric model (in the Kondo regime), switching 
on Ucf results in the increased itinerancy of the / elec¬ 
trons, however, above a certain value of Ucf it tends to 
localize / electrons again. The itinerancy of the conduc¬ 
tion electrons is gradually reduced as Ucf is increased. 
These results agree well with what has been obtained 
by quantum information analysis.— We also investigated 
what happens when the system is less than half-filled, 
that is the ground state is metallic. It has been revealed 
that in the mixed valence regime Ucf makes both the con- 
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duction and / electrons more correlated. These findings version of the DMRG methodi^^— 
agree qualitatively well with the properties of the infi¬ 
nite dimensional model, although, as expected for a one¬ 
dimensional model, there is no sharp Fermi edge. This ACKNOWLEDGMENTS 
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